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HIGHER ORDER APPROXIMATIONS OF 
RUNGE-KUTTA TYPE 

SUMMARY 


The fourth-order formula of Runge-Kutta has been used extensively to obtain 
approximate solutions of differential equations of first, second, and higher orders. This 
formula requires only four evaluations of the function involved in the differential equation 
while the fifth-order formula of Kutta-Nystrom requires six evaluations of the function. 
The sixth-order formula of Huta is not known as well, and has the disadvantage of re- 
quiring eight evaluations of the function. No formulas of higher order except those de- 
veloped in this investigation seem to exist in the literature. These formulas are "self- 
starting" and therefore extremely useful in combination with efficient "continuing" 
procedures. 

Formulas of higher order, including the seventh and eighth orders, are developed 
in this paper. Moreover, more efficient formulas requiring fewer evaluations are de- 
veloped for orders five and six. All possible formulas are displayed for those orders 
below five. 


FORMULATION OF THE PROBLEM 


1. General . The differential equations considered are of the form y' = f(x,y) , 
since systems of the first and higher orders are reducible to an application of this result. 
The function f is assumed to be analytic in a sufficiently large neighborhood of the initial 
point (x 0 , y 0 ) . 

Let f^ be defined by the equation 


i-i 

f i = f ( x o + a i h ’ y D + a i h £ b ij f j>’ 

j=0 

where f Q = f(x Q , y Q ) and i = i, 2, 


n 


Y = y 0 + h X c i f i 

i=0 


• • • 5 


n. Consider the finite series 



The a^, bj-:, Cj are parameters to be determined so that Y and the solution y(x 0 + h) of 
the given differential equation will agree to some desired degree of accuracy. 

It is immediate that a necessary and sufficient condition that the Taylor series 
for Y and the Taylor series for the solution y(x Q + h) of the given differential equation 
agree through terms in h™ is that 


-iMo 

1=0 


for k = 1, 2, m, where the symbol f^ is used for the k— derivative of the func- 
tion f(x Q + h, y(x 0 + h) ) . 

To specify a given formula, the coefficients may be detached as follows: 

Form ula 

a i a i b io 

a 2 a 2 b 20 + a 2 b 21 

# 

a n a n b no + a n b ni + • • • + a n b n,n-i 
Cq "t* c t + . . . 

This array specifies the formula 

n 

Y = y 0 + h Z c i f i » 

1=0 

where 

f o = f < x o-» yo> 

fj = f(x 0 + a. t h, y 0 + a t h b 10 f 0 ) 

f 2 = f(x 0 + a 2 h, y 0 + a 2 h(b 20 fo + b 21 fi) 

f n = f ( x o + a n h » y 0 + a n h ( b no f o + b ni f i + ••• + b n,n-i f n-i) ) 


2 



Such arrays will be used consistently to specify such formulas. 


n i-i 

Let 1 = Yj c i and 1 = Yj for 1 = 2 , • • • > n * Let i t + 2i 2 + . . . + ki k = k, 

i=0 j=0 J 

Q* = a., and 





i-i 

2*i S 
J=i 





(2) 


for k = 2, 3 m-i. Finally, let 


XX TP 

i = (k+ 1) J °i Q.-,. (3) 

i=i 

for m > 2 and k = 1, 2, m-1. These equations are sufficient conditions that the 

equations (i) be satisfied, giving a formula of order m, as is proved in a paper pre- 
sented by the author at the summer meeting of the American Mathematical Society in 
Boulder, Colorado, on August 29, 1963, An abstract of this paper appears in the 
American Mathematical Society Notices, Volume 10, Number 5, Issue No. 69, August, 
1963. The subscripts t and t r in the are introduced to distinguish the different 
products for fixed i and k. 

The equations ( 3) can now be written for each k = 1, 2, ..., m-1 in succession. 
Such equations for the first two values of k are given next. 


k = 1: 


1 = 2 X °i, a i, 

It— 1 


k = 2 : 


n 


Cj a? 


1 ' 3 j,-i ' J > “ 11 


1 = 6 


n h-i 

i -1 \i?i bj2jiaJl 

J2 -1 VJl— 1 j 
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In general, a simple induction establishes the facts that in each instance the 
indices for each by are distinct and the second index j is summed for i, 2, . , i-1, 
there are no other changes of indices, and each b-- is followed immediately by at least 
one aj . Hence, without confusion, b- a^ may be denoted by b K and the summation 
signs and indices may be omitted, as for the first two values of k below. 


k = 1: 

1 = 2 

ca 


k = 2: 

1 = 3 

ca 2 

1=6 cab 


The equations (3) are completely determined when the Q. are known. The rule 

k lt; 

for the inductive determination of the Q given in (2) is given below in the "summation 
convention" notation. 




(h \ 

n 2 abQ' 

It 1 j 


• • • 



(4) 


Let a be multiplied by 1, 2 ab be multiplied (and summed) by each term of Q 1 , . . . , 
kab be multiplied (and summed) by each term of Q^“^. Then the rule (4) states that a 
term of Qk is formed by forming a product consisting of if factors of the type mentioned 
first in the preceding sentence, i 2 factors of the type mentioned second in the preceding 
sentence, . . . , and i^ factors of the type mentioned last in the preceding sentence, where 
ij + 2 i 2 + . . . + k i^ = k. 

It is possible to state a rule so that the equations for a given k can be written 
without reference to the equations for smaller integers. Such a rule follows essentially 
from the following facts, which are easily proved: ( 1) is a product of factors, each 
factor of which consists of a sequence of an integer, a's, and b's, (2) in each factor of 
Qk, at least one a precedes each b and at least one a follows each b, (3) in each factor 
of Qk, the integer p immediately precedes a^b if and only if p - 1 is the number of a's 
that follow and are connected by summation to the given b, and (4) Qk contains exactly 
the number k of a's. However, such a rule is not of much advantage, if any, since all 
equations for k = 1, 2, . . . , m-1 are conditions to be satisfied to obtain a valid formula 
of order m. 


To further simplify the notation, let 
n-k 


J cj k) a. by = c! k+1) for j = 1, 2 n-k-i, 

i=j+i J 


(5) 
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where k = 0, 1, . .. , m-3 and, when convenient, c\° , c\ , c v . , ... will be denoted 
by c^, c^, c f . f , . . . respectively. Among the sufficient conditions determined later are 


1 = ( j + i) ( j + 2 ) . . . ( j + k + 1) 


n-k 

E 


i=i 




( 6 ) 


where j = i, 2, ... , m-k-1 and k = 0, i m-2. 

/ k) 

For m = n + i, it is seen that c x l ^ 0 from the definition and the fact that 
. n-k . 

1 = (n + i) I cj n a 1# Hence equations (5) determine the b^ in terms of the c . and 

(k) J 1 

the a^. Moreover, equations (6) determine the c\ in terms of the a.^ In this case, 

the algebraic system is reduced to the conditions not included in equations ( 6) , given the 
n free parameters a^. For example, when m = n + 1 = 4, there is one such condition 

i = 8 ca 2 b and three free parameters a ls a 2 , a 3 . 

A listing is given below of the for k = i, 2, 3, 4. Formation of the list is 

made easier by separating the products of order p + i of the type pabQ^ ^ from the 
other products of order p + i as is done below. 


k = i: 

a 

<NI 

II 

a 2 

k = 3: 

a 3 2a 2 b 

II 

a 4 2a 3 b 3a 2 b 2 


6a 2 bb 4a 2 (b) (b) 


2ab 

3ab 2 

4ab 3 

24abbb 


The term 4a 2 (b) (b) is formed by the rule from ( 


6abb 

8ab 2 b 12abb 2 


2ab) ( 2ab) and is written in summation 


notation 




i 2 -i 





b. 


L=i 


hh 



or simply 
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4a? £ b . . a . . 

12 V i t =i 1211 11 


The term 6a^bb is written 


6a a X b. . a. £ b a. 

■> fa V2 l ^ Vl >1 


in summation notation. 

For k = 4, it can be shown directly by an inductive process that the two conditions 
1 = 30 cabb and 1 = 40 c'a 2 b can be replaced by the one condition 7 = 120(cabb + c'a 2 b) 
and that the resulting fifteen conditions for k = i, 2, 3, 4 are necessary and sufficient 
for obtaining a valid formula for order 5. This was shown for i ^ k < 5 ( which includes 
the case i < k < 4) in a paper by Huta [1]. The necessary and sufficient conditions for 
1 < k < 6 have been derived by the author and will be given below, although the sufficiency 
of the other conditions are generally enough for our purposes. 


?r 

ii 

1 = 2ca 


<N| 

II 

1 = 3ca 2 

1 - 6c'a 

k = 3: 

1 = 4ca 3 

1 = 12c'a 2 

FT 

II 

1 = 5ca 4 

1 = 20c 'a 3 


1 = 10ca 3 b 

1 = 20ca 2 (b) ( b) 

k = 5: 

1 = 6ca 5 

1 = 30c'a 4 : 


1 = 12ca 4 b 

1 = 24ca 3 ( b) (b) 


13 = 720ca 2 (b) (bb) +360c’a 2 (b) (b) 2 = 45(ca 3 bb + c'a 3 b) 

1 = 360c ,M a 2 1 = 36ca 2 (b) (b 2 ) i = 40(ca 2 bb 2 + c'a 2 b 2 ) 

1 = 720c""a 1 = 60(ca 2 bbb + c’a 2 bb + c"a 2 b) 


ii i mi i mu iiib in mil 



k = 6: 


i = 42c'a 5 


1 = 35ca 2 b 4 


i - 
1 = 
i = 
5 = 
3i = 


19 = 


1 = 
1 = 
1 = 
1 = 
1 = 
19 = 


7ca 6 


1 = 210c"a 4 


14ca 5 b 1 = 21ca 4 b 2 1 = 28(ca 4 bb + c'a 4 b) 

28ca 4 (b) (b) 1 = 28ca 3 b 3 9 = 280(ca 3 b 2 b + ca 2 b 3 b) 

252( ca 3 bb 2 + c'a 3 b 2 ) 1 = 42ca 3 (b) (b 2 ) 1 = 105ca 2 b 2 b 2 

2520( ca 3 bbb + c'a 3 bb + c"a 3 b) 5 - 336ca s (b) (bb) + 168c'a 3 (b) (b) 
1680(ca 2 b 2 bb + ca 2 bb 2 b + c’a 2 b 2 b) 11 = 840(ca 2 bb 3 + c'a 2 b 3 ) 
80ca 2 (b) ( b 2 b) + 40ca 2 b 2 ( b) ( b) 1 = 56ca 3 ( b) ( b) ( b) 

16 8 ( ca 2 bbb 2 + c’a 2 bb 2 + c"a 2 b 2 ) 1 = 63ca 2 (b 2 ) (b 2 ) 

280(ca 2 bbbb + c'a 2 bbb + c"a 2 bb + c'"a 2 b) 

56[ca 2 (b) (bb 2 ) +ca 2 (b 2 )(bb) + c'a 2 (b) (b 2 ) ] 

56ca 2 (b) (b 3 ) 1 = 840c"’a 3 1 = 2520c’"'a 2 1 = 5040c"’"a 

1260 [2ca 2 (b) (bbb) +ca 2 (bb)(bb) + 2c’a 2 (b) (bb) +c”a 2 (b)(b)] 


This list could be extended, although the computation would be laborious. In- 
stead, the extension of the list giving sufficient conditions is relatively easy. For k = 5, 
the two lists agree with the exception that certain formulas in the list just given should 
be replaced with other formulas as shown next. 


13 = 720ca 2 (b) (bb) + 360c’a 2 (b) (b) 

1 = 72ca 2 (b) (bb) 

1 = 120c’a 2 (b) (b) 

2 = 45(ca 3 bb + c'a 3 b) 

1 = 36ca 3 bb 
1 = 60c’a s b 

1 = 40(ca 2 bb 2 + c'a 2 b 2 ) 

1 = 72ca 2 bb 2 
1 = 90c’a 2 b 2 

1 = 60(ca 2 bbb + c'a 2 bb + c"a 2 b) 

1 = 144ca 2 bbb 
1 = 180c'a 2 bb 


1 = 240c"a 2 b 
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The procedure is the same for k = 6 and is clear from our previous rules. 

If m = n + i and a =1, then it follows, for each i, that 
n ’ 


c i 


c i a i 


c. 

1 


( 7 ) 


where c^ = 0 is taken. In general , even though m and n + 1 are different integers , we 
may choose a n = 1 and assume the validity of equations (7) for each i. Since an obvious 
requirement is that m < n + 1, this will never lead to a contradiction per se. Under 
these circumstances, a great simplification occurs in the conditions for k = 1, 2, . .. , 
m-1. Most equations involving c^ can be omitted because of equivalent equations in- 
volving c^. Some exceptions are equations like 1 = ca and 1 = 20ca 2 (b) (b) . Further 
important simplifications can be made by requiring some of the following conditions to 
be satisfied for appropriate i. 


i-i 


i-l 

i-l 

Y. b.. a. 

= a. 

3 y b.. a? = a? 

6 z 

J=i 1] J 

i 

i=l 1J ] 1 

j=l 

i-1 

A b ij a j s 
i-i J 

= a? 

i 

• • • 



k=l 


jk a k “ a i 


However, this will be treated in detail in the discussion of formulas for a given order. 


GENERAL FORMULAS FOR ORDERS BELOW 4 


The formulas in these cases are easy to derive and are simply listed below. 

m = 0: Y = y 0 

m = 1: Y = y 0 + hf 0 

m = 2: Y = y 0 + h ^ [ ( 2a j - 1) f 0 + fj ] 
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m = 3: i. Y = y 0 + h [f 0 + ( 3 - 6c) f t + 6cf 2 ] 

h = f(x 0 +| h, y 0 +| hf 0 ) 

f 2 = f(x 0 +| h, y 0 + h^- [(4c - i)f 0 + fj) 

2. a^ a^ 

a 2 a 2 (t >20 + t> 2 i) 

(00 + 0! + C 2 ) 


where 

t»2i 


Cl 


a 2( a 2 ~ a l) 
a i( 2 - 3aj) 

2 - 3a 2 
(3 a i( a i - a 2 ) 


^20 ~ * “ ^21 


_ 2 - 3a< 

° 2 “ 6a 2 (a 2 - a t ) ’ 


Cq — i Cj - C2 


For the last formula, we have detached the coefficients as explained previously. This 
will be done for subsequent formulas. 


GENERAL FORMULAS FOR ORDER 4 


In this case, n = 3 is taken which is its minimum value. Denote 


- i-i ± 

b. = Yj b.. a. - — 

1 j=± 3 2 


"2 a i 


( 8 ) 


for i = 1, 2, 3, where b t = - ~r . Then three of the seven conditions may be written 


cab = 0, ca 2 b = 0, c'ab = 0. 
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As vectors in three space, c, ca, and c’ are orthogonal to the non-zero vector ab and 
hence are linearly dependent. Thus, for i = i, 2, 3 


A Cj + B c*a. + Cc! = 0 

1 xi i 

for some non-zero vector (A, B, C) . Multiplying and summing in turn by a^ and a? , 
we get 


i A+ i B+ i c = °> 


i A + 4 B + 12 C = °- 


These equations imply, for each i. 


c. - c.a. = c' and a, = i since Co + 0. 
l x l i 6 d 


With this information, it is not difficult to derive the complete solution for order 4. 


First, a 2 + i. Otherwise, c£ = 0 and Cj = c 2 a 2 b 21 = 0, but this contradicts the 
requirement that 1 = 24ci a t . Then there are three cases as discussed next. 

i. a t = a 2 . Then 1 = 6(c} + c 2 )a 2 and i = 12(cJ + c 2 ) af. Hence a 1 = a 2 = 
and c} + c 2 =-^ . Let c 2 = 2c. In this case, the formula is 

O 


1 

2 


1 

2 


1 

2 


1 

i2c 


[(6c - i) + 1] 


i [0 + (1 - 6c) + 6c] 

^ [1 + (4 - 12c) + 12c + 1] 


The case 1 = 6c is the famous Runge-Kutta Fourth Order Formula. 
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to |m* 




2. a t = i. 

Then c 2 = 


Then = 0, 6c 2 a 2 = i, 12c 2 a 2 = i and it follows that 
1,2 . 1 111 
T~T 2 ° 2= i’ sothat <H + c 3 =--c 2 a 2 


a 2 - 
Let f 3 = c. 


In this case, the formula is 


1 

1 

2 

1 


3. 


1 

g ( 3 + 1) 

[(12c - 3) + (-1) +4] 

•|[l + (l-6c) + 4+ 6c] 
b 

a lf a 2 , 1 are distinct. In this case, 


cl' = 


1 - 2a 9 


1 24aj ’ 


ci = 


12a 1 (a 1 - a 2 ) 


c 9 = 


1 - 2a, 
12a 2 (a 2 - a t ) 


and the formula is 


a 2 

1 


where 


a 2( t>20 + b 2l) 

( b 30 + b 31 + b 32) 

(C 0 + Cj + C 2 + C 3 ) , 


b 21 “ 


a 2 ~ a l 

2a t ( 1 - 2a 4 ) 


t>20 “ 


3a t - 4a* - a 2 
2a 4 ( 1 “ 2a 


b 31 “ 


( 1 - a 4 ) ( 2 - a t - 5a 2 + 4a 2 ) 
2a 1 (a 1 - a 2 ) ( 3 - 4a t - 4a 2 + 6aja 2 ) 


b 32 


( 1 - H) ( 1 ~ a 2 ) ( 1 - 2a t ) 
a 2 (a 2 - ^j) ( 3 — 4a t 4a 2 + 6a^a 2 ) 
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to |m* 



£>30 - * - t>3i - b 32 


2a 1 - 1 


12a 1 (a 1 - a 2 )( ai - i) 


12a 2 (a 2 - a^) (a 2 - 1) 


3 - 4a 4 - 4a 2 + 6a 4 a 2 
12( ai - 1) (a 2 - i) 


C 0 — i - Cj - c 2 - Cg 


FORMULAS FOR ORDER 5 


It will now be shown that the fifteen conditions for order 5 are inconsistent when 
n = 4. the minimum value. 


Six of the conditions may be written 


cab = ca 2 b = ca 3 b - c'ab = c"ab = (ea 2 b + c'a^b = 0 , 


where the notation (8) is used. In four space, the vectors c., c.a., c.a?, c!, c" 

^ J ' J J J J 

and /, c.a. 2 b.. + c'. a. for j = 1, 2, 3, 4 are orthogonal to the non-zero vector ab . 

.V, i 1 ij ] J 

i-j+1 ^ 

Since c\ = c^'= c 3 ' = 0 by convention, AC| + Bc! + CcN 0 implies A = B = C = 0. 

J J J 

Hence c., c'., and c'.' are independent vectors, spanning the three space orthogonal to 
1 J J 

the space spanned by the vector ab . 

By a straightforward calculation, we get 


c^ = a 4 Cj + (7 - 8a 4 ) ct -20(1- a 4 ) cV 


c.a? = a|c. + (6 - 8aJ)c. - 2(9 - 10a|)cV 
J J J J J 


Y, c^a?b i4 + c^a 4 = 2c' - 3c'.' 
i=j+l 1 13 J 1 J J 

for j = 1, 2, 3, 4. From the last of these equations 


'4<±4 U 43 Tt-jdj- ^o 3 
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which, since c 4 a 4 b 4 3 = c 3 and c 3 + 0, is equivalent to the condition a 3 + a 4 = 2. From 
the first two equations 

c 3( a 3 _ a d = ( ? - 8a 4 )c 3 , 
c 3 ( a I " a I) = (6 - 8a|) c 3 , 

from which it follows that 

( a 3 + a 4 ) ( 7 - 8a 4 ) c 3 = ( 6 - 8a| ) c 3 . 

Since c 3 ^ 0 and a 3 + a 4 = 2, it follows that a 4 = i and hence a 3 = 1 also. Since 

c. - c.a. = c T . , we get cl = 0, a contradiction. 

] J J J 3 

This contradiction establishes the assertion that all of the conditions cannot be 
satisfied when m ~ n + i = 5. 

One procedure at this point is to increase n from 4 to 5, which introduces 
additional parameters. This was done by Kutta, whose formula of the fifth order with 
six evaluations of the function contained an arithmetic error, later corrected by Ny strom. 

However, it is not necessary to increase n from 4 to 5 if we adopt a new view- 
point. Instead of requiring that the conditions be satisfied, let us require only that they 
be approximately satisfied. In other words, if the condition is written in the form F = 0, 
then the new requirement is to be that F = e for a sufficiently small e. In particular, 
for a given order m, the possibility will be examined that all conditions below order m 
are exactly satisfied, while some or all of order m are approximated. 

For m = 5, let a 4 = 1, c. - c.a. = cl for i = i , 2, 3, 4, and b- = 0 for i = 2, 3, 4. 

Let the conditions below order 5 be satisfied. Then i = 6c T a = 6cab = 3ca 2 - 3c t a 2 , 
so that Cj = 0 since i = 3ca 2 . Hence c\ = 0. The new requirements then reduce to the 
following. 


1 = 6c'a 

1 = 12c'a 2 

1 - e' = 20c'a 3 

1 = 2ca 

i = 24c"a 

1 - e" = 60c"a 



1 - e"’ = i20c 


where e T , e M , and e MT are to be taken small. Since c{ = 0, 3e f = 3 - 5a 2 - 5a 3 + iOa 2 a 3 

in 

c l a i 

and since a 2 = 2b 21 a 4 = 2 — , we deduce a 4 (2 - 5a 2 ) = 2a 2 e Tf + 2e tM (a 1 - a 2 ) . 
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To recapitulate, there are six parameters e', e", e"', a 1 , a 2 , a 3 and five requirements, 
namely e', e", e'" are to be taken small and 

3e* = 3 - 5(a 2 + a 3 ) + 10a 2 a 3 , 
a t (2 - 5a 2 ) = 2a 2 e" + 2e'"(a 1 - a 2 ) . 

When these requirements are met, a valid formula may be obtained. 

The formula just developed is given by the following schema. 


a i a i 

a 2 a 2(^20 + ^ 21 ) 

a 3 a 3(t>30 + t>3i + b 32 ) 

1 (b 40 + b 41 + b 42 + b 43 ) 

(C 0 + C 1 + c 2 + c 3 + c 4) 


where 


Z 


H+l 



a i b ij 


o (k+1) 

J 


for j = i, 2 


3-k and k = 0, 1, 2 ; b^ 0 = 1 - b^ 


- b. . for i = 2, 3, 4; and 
1,1-1 


120 a 4 


- 2e" + 2 - 5a, 

TT 

-2e" + 2 - 5a 4 


°i 120a 1 (a 1 - a 2 ) 

c 2 - 

120a 2 (a 2 - a 4 ) 


O 

II 

T 

1 - 2a 3 

1 - 2a 2 

c 2 “ 

i2a 2 (a 2 -a 3 ) 

12a 3 (a 3 - a 2 ) 

o 

II 

o 

c 2 = 

i - 2a 3 

1 - 2a 2 

12a 2 (1 - a 2 )(a 2 - a 3 ) 

12a 3 (i - a 3 ) (a 3 
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C 0 — i - C| C 2 Cg 


3e' = 3 - 5(a 2 + a 3 ) + 10a 2 a 3 
a 4 (2 - 5a 2 ) = 2a 2 e" + 2e"’(a 1 - a 2 ) 
e', e", e ,M sufficiently small. 

It is clear that the requirements can be satisfied and that the schema defines a four 
parameter set of valid formulas of the fifth order with only five evaluations. Each such 
formula does not qualify in the traditional sense to be called a fifth order formula, but 
does qualify as a fifth order formula from the new viewpoint and from a practical stand- 
point. Specific formulas for this and other cases are given in a later section. 

FORMULAS FOR ORDER 6 


In this case, let a 5 = 1, c 4 - cia 4 = c| for i = i, 2, 3, 4, 5, and. b 4 = 0 for 
i = 2, 3, 4, 5. Let the conditions below order 6 be satisfied. Then 1 = 120c ,M a = 
120c”ab = 60c"a 2 - GOc’j'a 2 , so that Cj' = 0 since i = 60c"a 2 . Similarly, c 4 = c 4 = 0. 
The new requirements then reduce to the following. 

i = 6c'a i = 12c'a 2 1 = 20c'a 3 i - e' = 30c'a 4 

1 = 24c”a 1 = 60c"a 2 i - e" = 120c n a 3 

1 = 120c M, a i - e ,M = 360c ,M a 2 
i - e"" = 720c""a 

i = 2ca i - e = I80c’a 2 bb i - d = 90c’a 2 b 2 , 


where d, e, e', e", e' M , e"" are to be taken small. 

Since cj = c 4 = 0, it follows that 

2e’ = 2 - 3(a 2 + a 3 + a 4 ) + 5(a 2 a 3 + a 2 a 4 + a 3 a 4 ) - 10a 2 a 3 a 4 
e" = i — 2(a 2 + a 3 ) 5a 2 a 3 . 


( 9 ) 
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c , 1 '"a 1 


, it follows that 


III!!! 


Since a 2 =2b 21 a 4 


' 2' a 2 


a 4 (l - 3a 2 ) = e'"a 2 - e""(a 2 - a 4 ) . 


( 10 ) 


Laborious calculation leads to the 


1 - 180c'a bb = 


. 2 , , ( 1 - a 4 ) ( i 

I q = 


„ „ ( i - a 4 ) ( 2 

1 - 90c'a 2 b 2 = 4L1 — 


following results. 

- 5a 2 ) + e""( 3 - 5a 2 - 5a 4 + 10a 2 a 4 ) 

2(2- 5a 2 ) 

- 10a 2 + 15a 4 a 2 ) + 2e"'(3 - 5a 2 - 5a 4 + 10a 2 a 4 ) 

4(2 - 5a 2 ) 


From the equations (9) , we get 


1 - a 4 


2e' + e" 

3 — 5a 2 — 5a 3 10a 2 a 3 


which implies 1 - a 4 , i - 180c’a 2 bb, and 1 - 90c'a 2 b 2 are small provided e', e" , e'", 
e"" are small. Hence equations 1 - e = 180c'a 2 bb and 1 - d = 90c'a 2 b 2 will be satis- 
fied in the presence of the other requirements. 

The formula just developed is given by the following schema. 


a l a i 

a 2 a 2( t*20 + t>2l) 

a 3 a 3(t*30 + + b 32 ) 

a 4 a 4 (b 40 + b 41 + b 42 + b 43 ) 

i ( b 50 + b 51 + b 52 + b 53 + b 54 ) 

(c 0 + c i + c 2 + c 3 + c 4 + c 5 ) 
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where 


Z 


i=j+i 




c 


(k+i) 

3 


for j = i, 2, . . . , 4 - k and k = 0, 1, 2, 3; bj. = 1 - b t - . . . - b. . for i = 2, 3, 
4, 5; and ’ 1 * 


C HH 

1 


1 - e fm 
720a 1 


C MT 

1 


1 - e™ - 3a 2 
30Oa 4 (a 4 - a 2 ) 


1 - e"’ - 3a H 
360a 2 (a 2 - a t ) 


c\' = 0 


2 - 5a, 

c" = — — — 

2 120a 2 ( a 2 - a 3 ) 


c" 


3 


2 - 5a, 

120a 3 ( a 3 - a 2 ) 


i i 3 ~ 5a^ 5a 4 + 10a 3 a 4 ^ 3 ~ 5a 2 — ^^4 ^ 10a 2 a 4 i 3“5a 2 —5a 3 “i"i0a 2 a 3 

1 2 60a 2 (a 2 - a 3 ) (a 2 - a 4 ) 3 60a 3 (a 3 - a 2 ) (a 3 - a 4 ) * 60a 4 (a 4 -a 2 )(a 4 -a 3 ) 

3 5a 3 * 5a 4 + iOa 3 a 4 3 5a 2 — 5a 4 + 1 0aoa 4 

1 2 60a 2 <i - a 2 )(a 2 - a 3 )(a 2 - a 4 ) 3 60a 3 ( 1 - a 3 ) (a 3 - a 2 ) (a 3 - a 4 ) 

3 - 5a 2 - 5a 3 + 10a 2 a 3 
° 4 ~ 60a 4 ( 1 - a 4 ) (a 4 - a 2 ) (a 4 - a 3 ) 

i 

C 5 ~ £ “ c 2 a 2 “ c 3 a 3 “ c 4 a 4 c 0 ~ 1 ~ c 2 “ c 3 “ c 4 “ c 5 


2e T = 2 - 3(a 2 + a 3 + a 4 ) + 5(a 2 a 3 + a 2 a 4 + a 3 a 4 ) - 10a 2 a 3 a 4 


e” = i - 2(a 2 + a 3 ) + 5a 2 a 3 

a 4 ( 1 - 3a 2 ) = e M, a 2 - e'"'(a 2 - a t ) 

e' , e", e M, } e M " sufficiently small 
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It is clear that the requirements can be satisfied and that the schema defines a five 
parameter set of valid formulas of the sixth order with only six evaluations (in the 
sense of the new viewpoint) . 


EFFICIENT FORMULAS THROUGH ORDER 8 


Using the principles already explained, we have derived other formulas. To 
avoid going into details and to conserve space, we will simply list a spectrum of such 
formulas, including formulas of the seventh and eighth orders. 

Formula ( 4-4) 


1 1 
100 100 

( -4278 + 4425) 

—7- (524746 - 532125 + 16170) 
8791 v ' 

(- 179124 + 200000 + 40425 + 8791) 
70092 v 


3 

5 

1 


Formula ( 5-5) 


1 

9000 

_3_ 

10 

3 

4 

1 


1134 


1 

9000 

-i- (_ 4047 + 4050) 

■| (20241 - 20250 + 15) 

~7 (- 931041 + 931500 - 490 + 112) 

81 

(105 + 0 + 500 + 448 + 81) 


( 11 ) 


(12) 
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1 

300 

i 

5 

3 

5 

14 

15 

1 


Formula (6-6) 


1 

300 

4 (- 29 + 30) 

5 

~ (323 - 330 + 10) 

5 

— J— ( -510104 + 521640 - 12705 + 1925) 

( - 417923 + 427350 - 10605 + 1309 - 54) 
(198 + 0 + 1225 + 1540 + 810 - 77) 

3696 


Formula ( 7-7) 


1 

192 

i 

6 

1 

2 

1 

5 

6 

1 


1 

192 

| (- 15 + 16) 

~~ (4867 - 5072 + 298) 

1 86 

~7 (- 19995 + 20896 - 1025 + 155) 

« 3 1 

— (- 469805 + 490960 - 22736 + 5580 + 186) 

50z2 

— 7— (914314 - 955136 + 47983 - 6510 - 558 + 2511) 
2604 

(14 + 0 + 81 + 110 + 0 + 81 + 14) 

oUU 


(13) 


(14) 
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<N|0} | CO •'H | (N -rH|CD 00|C5 tH | 05 LQ | CO ^ | (Nl CSJ j C5 tH j CO -tH | CM <M | CO 


Formula (7-9) 


1 

2140320 


2 

9 


li5< 1 + 3 > 


| (1 + 0 + 3) 


216 


(23 + 0 + 21 - 8 ) 


729 


(- 4136 + 0 - 13584 + 5264 + 13104) 


(15) 


(105131 + 0 + 302016 - 107744 - 284256 + 1701) 

151632 

— — (-775229 + 0-2770950 + 1735136 + 2547216 + 81891 + 328536) 

1375920 v 


251888 


( 23569 + 0 - 122304 - 20384 + 695520 - 99873 - 466560 + 241920) 


( 110201 + 0 + 0 + 767936 + 635040 - 59049 - 59049+635040 + 110201) 


Formula (8-10) 


27 


Ii< 1 + 3) 


12 (1 + ° + 3) 


( 16 ) 


^(1 + 0 + 0 + 3 ) 


(13 + 0 - 27 + 42 + 8) 
54 
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1 

6 

i 

5 

6 

1 


( 389 + 0 - 54 + 966 - 824 + 243) 

(- 231 + 0 + 81 - 1164 + 656 - 122 + 800) 

(- 127 + 0 + 18 - 678 + 456 - 9 + 576 + 4) 

288 v 

( 1481 + 0 - 81 + 7104 - 3376 + 72 - 5040 - 60 + 720) 
(41 + 0 + 0 + 27 + 272 + 27 + 216 + 0 + 216 + 41) 


(16) 
( Cont'd) 


i 

9 

1 

6 

1 

4 

J_ 

10 

1 

6 

1 

2 

2 

3 

i 

3 

5 

6 


Formula (8-12) 

1 

9’ 

■h < 1 + 3 > 

1i< 1 + 0 + 3 > 

^<29 + 0 + 3 3 -i2) 

(33 +0 + 0 + 4+ 125) ( 17) v 

(- 21 + 0 + 0 + 76 + 125 - 162) 

36 

(- 30 + 0 + 0 - 32 + 125 + 0 + 99) 

Zt:0 

( 1175 + 0 + 0 - 3456 - 6250 + 8424 + 242 - 27) 

ijZrc 

■— (293 + 0 + 0 - 852 - 1375 + 1836 - 118 + 162 + 324) 
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5 1 

6 1620 ( 1303 + 0 + 0 " 4260 “ 6875 + 9990 + 1030 + 0 + 0 + 162) 

1 4428 ( ~ 8595 + 0 + 0 + 30720 + 48750 - 66096 + 378 - 729 - 1 944 - 1 296 + 3240) 

(41 + 0 + 0 + 0 + 0 + 216 + 272 + 27 + 27 + 36 + 180 + 41) ( 17) 

( Cont’d) 

Formulas (11) , ( 15) , and (17) have the orders indicated in the traditional sense while 
formulas ( 12) , ( 13) , ( 14) , and ( 16) have the orders indicated from the new viewpoint 
and only for certain values of h. 
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